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THE EFFECT OF ERRORS IN THE ASSIGNMENT 
OF THE TRANSMISSION FUNCTIONS ON THE ACCURACY 
OF THE THERMAL SOUNDING OF THE ATMOSPHERE 

Yu. M. Timofeyev, 

USSR Academy of Sciences 

In accordance with a program of Soviet-American cooperation /l* 

in the improvement of the methods of thermal sounding of the 
atmosphere, the section of atmospheric physics of the NIIF of 
the LGU performed an investigation of the effect of errors in 
the assignment of transmission functions on the accuracy of 
the indirect method. These studies included calculations of 
the "effective background noise" of measurement and numerical 
experiments on solving the problem in reverse by replacement. 

The corresponding calculations and numerical experiments were 
performed for conditions which functionally modeled the spe- 
cialized satellite-borne radiometers of the USSR and USA in 
space used for thermal sounding of the atmosphere. The 
spectral distributions of the measurement channels of thermal 
radiation are presented below. 

USSR Radiometer 


Channel 
Number N 

I 

2 

3 

4 . 

5 

6 

Central 

Frequency (cm -1 ) 

666,5 

676,0 

692,5 

699,0 

736,0 

746,0 


* Numbers in the margin refer to the pagination in the foreign text. 
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USA Radiometer 


Channel 
Number N 

I 

2 

3 

/ 4 ' 

5 

6 

7 

Central 

frequency (cm -1 ) 

669,0 

679,0 

690,0 

702 ,P 

719,0 

732,0 

751,0 


In order to obtain a model of the errors in the trans- 
mission function of the atmosphere (and consequently in the 
root of the corresponding integral equation [1,2]), the 
transmission functions were obtained by means of a program /2 

of linear calculations for the following model atmospheres: 

1. The standard model of the atmosphere with a carbon 
dioxide content c = 330 ppm ("basic" variant). 

2. A standard model of the atmosphere with c = 363 ppm 
(the "content" variant) . 

3. A standard model of the atmosphere (c = 330 ppm) 
with a displacement of the function of the measurement channel 
of the apparatus by 1 cm" 1 in the axis of frequency (the 
"displacement" variant) . 

4. A standard model (c = 330 ppm) with an increase in 

t he half-width of the function of the channel of the apparatus 
by 10% (the "stretched" variant) . 

5. A standard model (c = 330 ppm) including "exponential" 
type aerosol absorption and an integral optic depth T a (0,°°) = 

0.1 (the "aerosol 1" variant). 

6. A standard model (c = 330 ppm) taking into account 
"layered" aerosol absorption and an integral optic depth 

= 0.1 (the "aerosol 2" variant). 

7. The sanu> as variant 5, but with T a (0,°°) = 0.2. 

8. The same as variant 6, but with r a (0,») = 0.2. 



The calculations thus performed made it possible to 
obtain the variation in the transmission function PAv and 
in the roots of the integral equation, relative to the "basic" 
variant, for the USSR and USA radiometers, caused by variations 
in the carbon dioxide abundance, errors in assigning the 
spectral characteristics of the satellite apparatus, and the 
presence of additional atmospheric absorption (in the aerosol 
frequencies) . 

Table 1 presents the maximum values of the variation of 
the function PAv for various variants of calculation with 
respect to the "basic" variant. The value of maximum variation, 
as the data presented here imply, does not exceed 0.1; that is, 

10% absolute transmission. The maximum variation, APAv, is 
observed in the first channel of the USA radiometer (at the 
center of the band of absorption) with a displacement of the 
apparatus function (AP = 9.9%) and in the "transparent" channels 
of the USSR and USA radiometers at maximum aerosol optical depth /3 
Ta(0*°°) = 0*2 (AP = 6.3 and 8.5%, respectively). 

Figure 1 illustrates the weighting functions [1] for six 
spectral intervals of the USSR radiometer as well as their 
variations caused by changes in the CO 2 concentration (fig. lb), 
"stretching" (fig. lc) , and "displacement" (fig. Id) of the 
apparatus function. It follows from figure 1 that although 
the vertical changes in the variation of the weighting functions 
are similar for all the mechanisms considered "disturbing" 
the transmission function of the atmosphere, the amplitudes 
of the variation vary considerably. 

Table 2 presents the values of the emitted radiation for /4 

the channels of the USSR and USA radiometers in the "basic" 
variant and the "effective noise" of measurement, AJ, repre- 
senting the variation in emitted radiation caused by the 
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TABLE 1 

MAXIMUM VARIATION OF THE TRANSMISSION FUNCTION 


a) USSR Radiometer 


■ — ■ ■ — 
Channel N 

Disturbance 

I 

2 

3 

4 

5 

6 

CO 2 

Variation 

0,036 

0,040 

0,041 

0,035 

0,011 

0,039 

Displacement 
4>(v-v' ) 

0,026 

0,007 

0,028 

0,026 

0,012 

0,019 

Stretch 

4>(v-v') 

0,008 

0,004 

0,003 

0,007 

0,006 

0,006 

Aerosol 1 
T a “ o.l 

0 

0 

0 

0,001 

0,017 

0,033 

Aerosol 1 

T a ■ 0.2 

0 

0 

0 

0,001 

0,033 

0,063 


b) USA Radiometer 


Channel N 






1 


Disturbance 

I 

2 

3 

4 

5 

6 i 

7 

CO 2 

Variation 

0,027 

0,041 

0,040 

• 

0,039 

0,038 


0,039 

Displacement 

4>(v-v') 

0,099 

0,012 

0,026 

0,029 

0,022 

0,006 

0,031 

Stretch 
4>(v— v * ) 

0,018 

0,003 

0,006 

0,006 , 

0,003 

0.006 

0,006 

Aerosol 1 
T a * 0.1 

0 

0 

0 

QH 

0,007 

0,017 

0,045 

Aerosol 1 
T a - 0.2 

0 

0 

0 

0,001 

0,014 

0,033 

0,085 
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TABLE 2. EMITTED RADIATION AND "EFFECTIVE NOISE" 
a) USSR Radiometer (erg/cm 2 at average cm" 1 


Channel 
Number N 

Disturbance 

I 

i* 

2 

3 

4 

5 

6 

Basic 

Variant 

52,14 

47,96 

44,60 

50,77 

75,87 

88,32 

Content 

0,70 

0,62 

0,15 

1,22 

2,83 

2,12 

"Stretch" 

0,039 

0,0041 

0,12 

0,35 

0,10 

0,13 

Displacement 

0,73 

0,23 

0,24 

0,88 

0,71 

1,14 

Aerosol 1 
la = 0.1 

0,0020 

0,0002 

0,0045 

0,095 

1,01 

1,72 

Aerosol 1 
Ta = 0.2 

0,0041 

0,0005 

0,0090 

0,18 

1,96 

3,34 

Aerosol 2 
la * 0.1 

0,0095 

0,0037 

0,022 

0,19 

1,35 

2,13 

Aerosol 2 
T a “ 0.2 

0,0023 

0,0094 

0,068 

0,43 

2,75 

4,27 


b) USA Radiometer (erg/cm 2 at average cm" 1 ) 


Channel 
J^umber N 

Disturbance' 

I 

2 

3 

r 

4 

5 

6 

7 

Basic 

Variant 

56,04 

46,37 

44,79 

51,57 

66,16 

77,73 

91,72 

Content 

0,98 

0,54 

0,059 

1,45 

2,59 

2,93 


"Stretch" 

0,47 

0,11 

0,17 

0,29 

0,15 

0,27 

0,3G 

Displacement 


0,18 

0,30 

1,20 

T 09 

0,41 

1,56 

Aerosol 1 

!j = 0.1 

0,00008 

0,0003 

0,0910 

0,091 

0,51 

1,05 

1,92 

Aerosol 1 
Ta ** 0.2 

0,00017 

fm 

0,0081 

0,17 

1,02 

2,94 

o. n/ i 

Of 

Aerosol 2 
T a • 0.1 


0,0041 

0,19 

0,20 

0,78 

1,41 

2,35 

Aerosol 2 
Ta “ 0.2 

0,0048 

0,011 

0,060 

0,47 

1,64 

2,86 

4,70 
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various variations in the transmission function (the variation 
of the radiation was calculated according to an average 
temperature profile proposed by the American side together 
with a corresponding covariant temperature matrix Kxx) • It 
follows from the data of table 2 that in the majority of cases 
the "effective noise" of measurement exceeds the level of 
average square accidental error oj of contemporary satellite 
apparatus and in particular exceeds the value taken for 
numerical experiments on the solution of the reverse problem 
(see below) oj = 0.15 erg/c at avg. cm” 1 . The exceptions are 
the values of the "effective noise" in the optically dense 
channels of measurement with aerosol disturbances (this is 
caused by the linear nature of the vertical changes in optical 
density i a (0,2)), and in the "stretched" apparatus functions 
<t>(v-v') for channels N = 1-3 in the USSR radiometer and N = 2 
in the USA radiometer. We should draw attention to the fact 
that the "stretched" 4>(v-v') value of AJ for the first two 
channels is substantially larger for the USA radiometer than 
for the USSR radiometer. This circumstance is a result of the 
spectral character of the first two channels of the USA radiom- 
eter than for the USSR radiometer. This circumstance is a 
result of the spectral character of the first two channels of 
the USA radiometer, in particular the noticeably smaller half- 
width <f>(v-v') compared to the USSR radiometer. From this, it 
also follows that a high accuracy of the assigned spectral 
character is necessary to a radiometer, especially in the area /6 
of rapidly changing spectral and optical properties of the 
terrestrial atmosphere. We should note, finally, that in the 
remaining cases the values of the "effective noise" are similar 
for the two radiometer types. 

Thus, the calculations of the "effective noise" of measure- 
ment in the analyzed variants of the transmission function 
show that these variations in the majority of cases (with the 
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possible exception of the "stretch" variant) appear to have a 
noticeable effect on the accuracy of the indirect determination 
of the temperature profile. 

During the next stage of the investigation, numerical 
experiments were performed on generating the vertical tempera- 
ture profiles by a replacement arrangement for 30 temperature 
profiles specifically chosen by the American side. For each 
temperature profile the emitted radiation JAv£ was calculated 
for the transmission function corresponding to the "basic" 
variant, the accidental errors Ej_ (distributed normally with 
a mean square deviation Oj = 0.15 erg/cm 1 2 at avg . cm -1 ) were 
added to the value JAv^, and then the quantity JAvi + ei was 
used to generate a vertical temperature profile. In this 
solution of the reverse problem, the transmission function 
corresponding to both the "basic" variants and to the other 
variants was used: the "content," "displacement," "stretch," 

and so on. 

The generation of the temperature profile took place by 
means of a method of statistical regularization [1,2] using 
the correlation matrix Ktt calculated by the American team 
from a global ensemble of 210 profiles. It is important to 
emphasize that the generation of each temperature profile was 
performed for a large set of accidental errors (50 examples), 
and that the average characteristics of the accuracy of the 
indirect method were subsequently analyzed. 

Let us enumerate the basic results and conclusions obtained 
in this stage: 

1. The errors in generating temperature varied widely /7 

for different profiles and different levels of the atmosphere. 

As an example, figure 2 presents the dependence of the mean 
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square error in the generated temperature, o T , in the 1000-10 mb 


levels, on the distance o K = .• »J 


w • I 


nV 


of the temperature from the average temperature profile. The 
figure illustrates the strong variability of o-j (from 1 to 4°K) 
and the weak correlation between oj and °K* 

2. A comparison of the accuracy of generating tempera- 
tures by means of the USSR and USA radiometers revealed a 
small advantage in the USA radiometer, possibly caused by an 
insufficient number of channels in the Soviet device for error 
calculations (oj = 0.15 erg/cm 2 at avg . cm -1 ). The mean square 
error oj (for 30 profiles, 28 levels in the range from 1000 to 
10 mb pressure and 50 instances of accidental error) of the 
Soviet (six-channel) radiometer consists of 2.20°K, and for 
the USA (seven-channel) radiometer, 2.05°K, at exactly assigned 
roots of the integral equation. This relatively low level of 
accuracy in the indirect method is explained by two basic 
causes : 


a) The ensemble of 30 temperature profiles used 
in the numerical experiments was chosen from an initial 
global ensemble of temperature profiles. According 1 7, the 
covariant matrix K TT constructed from the initial global 
ensemble (and used in the method of statistical regularization) 
describes the global structure of the thermal state of the 
atmosphere. In this connection, the a priori mean square 
temperature deviations at many atmo. heric levels exceeds 
10°K. Such a large a priori indeterminacy in the knowledge 

of the thermal structure at an essential point determines the 
accuracy of the indirect method. 

b) It is possible that six-channel and even seven- /8 

channel systems of measurement are not optimal for the level 

of accidental error (oj = 0.15 erg/cm 2 at avg. cm -1 ) used in 
the numerical experiments and for global observations. 
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3. The increase of the error? of generation created 
by errors in the assignment of roots to the integral equation 
is extremely changeable. These increases vary from profile 
to profile, are different at different levels of the atmos- 
phere, and so on. With an error in the C0 2 concentration of 
30% the mean square error increases in the USA radiometer by 
0.68°K and in the USSR radiometer by 0.52°K. It is interesting 
to note at this point that the effectiveness of both sounding 
systems is practically the same, which is easily explained by 
the fact that t total noise (the accidental error plus the 
"effective noise") has already reached the point where six- 
and seven-channel devices are practically equivalent. The 
maximum increase in error of generation AT m is observed close 
to the earth's surface (AT m = 3°K) and the region of the 
tropopause ( AT m = 1-2°K) . 

"Stretching" of the apparatus function of the radiometer 
(increasing the half-width by 10%) does not introduce any 
essential change in the error of generation. Thus when the 
mean square errors increase by 0.03-0.05°K, the maximum 
increase in errors reaches , v0.3 o K. 

A change in the central frequency of the channel (the 
systematic "displacement" of the apparatus function of the 
radiometer by 1 cm- 1 ) introduces an increase in the mean 
square error of generation o T of 0.21°K (the USSR radiometer) 
and of 0.27°K (the USA radiometer). The maximum increase for 
the USSR radiometer consists of 1.63°K and for the USA radiom- 
eter, of 5.18°K, and it is observed in the upper levels of 
the atmosphere. The substantial increase in the error of 
generation T(p) in the upper atmospheric layers for the USA 
radiometer is caused by the peculiarity of the spectral char- 
acteristics of its first channel (see below) . 
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The effect of the aerosol absorption is determined above /9 

all by the value of the integral optical depth T a (0,°°). When 
the value x a (0,®) = 0.1, the increase in the mean square error 
is 0.16-0.18°K, and the maximum error (in the lower layers of 
the atmosphere) is 'v-3°K. The corresponding value at T a (0,°°) = 

0.2 are equal to 0.56-0.63°K and ^6°K. Table 3 presents the 
values of the increase in the mean square error (Ao?) and the 
maximum error (AT m ) in the generated temperatures for all the 
considered variants of "disturbed" transmission function. 


TABLE 3. INCREASE OF MEAN-SQUARE ERROR ( A o T ) AND MAXIMUM 
ERROR ( AT m ) OF GENERATED TEMPERATURE (IN °K) 


It- 

Ilodol 

Variant 

USSR Radiometer 

USA Radiometer 

aSt 

A i m 

AQr 

A i **• 

I 

"Content" 

I 0.52 

2.53 

0.68 

2.61 

2 ! 

1 j 

"Displacement" j 

0.21 

1.63 

0.27 

5.18 

3 

"Stretched" 

! 0.03 

0.33 

0.05 

0.31 

4 

Aerosol 1 j 

0.18 

2.64 

0.17 

2.80 

5 

Aerosol 2 j 

0,7.6 

2.46 

0.16 

2.62 

6 

Aerosol 3 

0,61 

6.53 

0,63 

6,75 

7 

Aerosol 4 

0.56 

6. II 

0.58 

6.32 

• 


Thus, all of the "disturbances" to the transmission function 
(except the "stretched") considered above substantially influ- 
ence that accuracy of the thermal sounding of the atmosphere. 
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On the basis of the results obtained, let us formulate 
the accuracy requirements of the assigned transmission function. 
In case we analyze the maximum increase of the error in gener- 
ated temperature AT m , we find that in order to ensure AT m does 
not exceed 1°K, it is necessary that the maximum variation in 
the transmission function does not exceed 'vO.Ol. An analysis 
of the relationship between the increase in the mean square 
error of the generated temperature Aox and the mean square 
errors of the assigned transmission function o p makes it pos- 
ible to conclude that, for o p < 0.02, Ao? usually does not 
exceed 0.5°K, at which point the maximum errors AP m can be 
0.03-0.04. As a qualitative average number, it appears 
advisable to infer an accuracy of assigning PAv to 2% abso- 
lute transmission. Given this level of error, the increase 
in the mean square, generated errors do not exceed an average 
of 0.5°K (or 25% of the initial value of o T ) . 


/10 
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LIST OF FIGURES 


Figure 1 


Figure 2 


The weight functions (upper left, measured in °K) 
corresponding to the measurement channels of a 
Soviet radiometer, plotted against pressure P in 
millibars; the deviations (measured in A°K) caused 
by the variants in C0 2 content (upper right) , 
increasing half-width (lower left) , and displace- 
ment (lower right), plotted against pressure. 


Plot of the mean square errors of generation ox 
(y-axis, in °K) versus the values of initial 
deviation (x-axis, in °K) for 30 temperature 
profiles. 
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